Thick planar domain wall: its thin wall limit and dynamics 
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Abstract 

^f} ' We consider a planar gravitating thick domain wall of the Xcj) 4 theory as a spacetime 

with finite thickness glued to two vacuum spacetimes on each side of it. Darmois junction 
conditions written on the boundaries of the thick wall with the embedding spacetimes repro- 
duce the Israel junction condition across the wall in the limit of infinitesimal thickness. The 
thick planar domain wall located at a fixed position is then transformed to a new coordinate 
system in which its dynamics can be formulated. It is shown that the wall's core expands as 
if it were a thin wall. The thickness in the new coordinates is not constant anymore and its 
time dependence is given. 
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1 Introduction 



Domain walls are solutions to the coupled Einstein-scalar field equations with a potential having 
a spontaneously broken discrete symmetry and a discrete set of degenerate minima. In the 
simplest case of two minima, a domain wall having a non-vanishing energy density appears in 
the separation layer, with the scalar field interpolating between these two values. 
Domain walls in the cosmological context have a long history pQ. It was realized very early 
that the formation of domain walls with a typical energy scale of > lMeV must be ruled out 
0, because a network of such objects would dominate the energy of the universe, violating the 
observed isotropy and homogeneity. Domain walls were reconsidered in a possible late time 
phase transition scenario at the scale of < WOMeV. Such walls were supposed to be thick 
because of the low temperature of the phase transition [Sj. The suggestion that Planck size 
topological defects could be regarded as triggers of inflation, revived the discussion of thin and 
thick domain walls [I]. The realization of our universe as a (3 + l)-dimensional domain wall 
immersed in a higher dimensional spacetime has led to the recent numerous studies pj. 
First attempts to investigate the gravitational properties of domain walls were based on the so 
called thin wall approximation. In this approach one forgets about the underlying field theory 
and simply treats the domain wall as a zero thickness (2+l)-dimensional timelike hypersurface 
embedded in a four-dimensional spacetime. The Darmois-Israel thin wall formalism [£j is then 
used to match the solutions of the Einstein equations on both sides of the wall in the embedding 
spacetime across the thin wall. However, such spacetimes have delta function-like distributional 
curvature and energy-momentum tensor supported on the hypersurface. Using this procedure, 
the first vacuum solutions for a spacetime containing an infinitely thin planar domain wall was 
found by Vilenkin j7., and Ipser and Sikivie [Sj. The very interesting feature of such domain 
walls is that they are not static, but have a de Sitter-like expansion in the wall's plane. External 
observers experience a repulsion from the wall, and there is an event horizon at finite proper 
distance from the wall's core. These results were initially obtained within the framework of 
the Darmois-Israel thin wall formalism which has been used as an approximative description 
of real wall of finite thickness. Typically, a self gravitating domain wall has two length scales, 
its thickness w and its distance to the event horizon which can be compared to w. Since these 
lengths are expressed in terms of the coupling constants of the theory, thin walls turn out to be 
an artificial construction in terms of these underlying parameters as mentioned in (Hj. 
The first exact dynamical solution to thick planar domain walls was obtained by Goetz ^U] 
and later a static solution was recovered with the price of sacrificing reflection symmetry 
Silveira studied the dynamics of a spherical thick domain wall by appropriately defining an 
average radius < R >, and then used the well-known plane wall scalar field solution as the first 
approximation to derive a formula relating < R > , < R > , and < R > as the equation of motion 
for the thick wall. Widrow ^3] used the Einstein-scalar equations for a static thick domain wall 
with planar symmetry. He then took the zero-thickness limit of his solution and showed that 
the orthogonal components of the energy-momentum tensor vanish in that limit. Garfinkle and 
Gregory Jl] presented a modification of the Israel thin shell equations to treat the evolution 
of thick domain walls in vacuum. They used an expansion of the coupled Einstein-scalar field 
equations describing the thick gravitating wall in powers of the thickness of the domain wall 
around the well-known solution of the hyperbolic tangent kink for a A(/> 4 wall and concluded 
that the effect of the thickness at the first approximation was to effectively reduce the energy 
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density of the wall compared to the thin case, leading to a faster collapse of a spherical wall 
as well as a slower expansion of a planer wall in vacuum. Barrabes et al. ^5] applied the 
expansion in the wall action and integrated it out perpendicular to the wall to show that the 
effective action for a thick domain wall in vacuum, apart from the usual Nambu term, consists 
of a contribution proportional to the induced Ricci curvature scalar. Arodz et al. ,10 applied a 
cubic polynomial ansatz for the scalar field inside the wall while the field gets the vacuum values 
outside the wall. Then the field equations and boundary conditions were used to determine the 
coefficients of polynomial. They obtained a Nambu-Goto type for the core of the wall together 
with a nonlinear equation for its thickness. Within the context of a fully nonlinear treatment of a 
scalar field coupled to gravity, Bonjour, Charmousis and Gregory (BCG) found an approximate 
but analytic description of the spacetime of a thick planar domain wall of the Xc/) 4 model by 
examining the field equations perturbed in a parameter characterizing the gravitation interaction 
of the scalar field ^7] , but they did not analyze the motion of the thick wall. The thin wall limit 
of Goetz's solution was studied in JH] and it was shown that this solution has a well-defined 
limit. 

A completely different method based on the gluing of a thick wall, considered as a regular 
manifold, to two different manifolds on each side of it was developed by Khakshournia-Mansouri 
(KM) in [12] and recently in a more generalized context in [20,. In this paper, we will first 
investigate the thin wall limit of the thick BCG-like domain wall adapted to the KM formalism. 
We then introduce a coordinate transformation from the BCG solution to the one in which one 
can observe evolution of the thick wall. 

The organization of this paper is as follows. In section 2 we give a brief introduction to BCG 
thick wall solution and summarize all the useful equations we need in this paper. Section 3 is 
devoted to the thick wall formalism developed in |19j and its application to the thick planar 
domain wall. The thin wall limit is then established in section 4. The dynamics of this thick 
domain wall in the new coordinates is studied in section 5. We then summarize our results in 
section 6. 



2 Planar domain wall solution 

Domain walls as regions of varying scalar field separating two vacua with different field values 
are usually described by the following matter Lagrangian: 

L = V^W-y^), (1) 

where (ft is a real scalar field and V(4>) is a symmetry breaking potential which we take it to 
be V(4>) = \{4> 2 — i] 2 ) 2 , where A is a coupling constant and r\ is the symmetry breaking scale. 
Looking for a static solution of the equation of motion derived from this Lagrangian in flat 
space-time, one gets 

A = tanh— , (2) 

w 

where X = w = an d z is the coordinate normal to the wall. This particular solution 
represents an infinite planar domain wall centered at z = 0. From the stress-energy tensor of 
the wall, one can easily observe that the wall energy density is peaked around z = and falls 
down effectively at z = w. We may therefore call w, being a characteristic length scale of the 
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domain wall, the effective thickness of the wall. 

We now look at the planar gravitating domain wall solutions. The line element of a plane 
symmetric spacetime may be written in the general form 



ds" 



A 2 {z)dt 2 - B 2 (z,t)(dx 2 + dy 2 



dz l 



(3) 



which displays reflection symmetry around the wall's core assuming to be located at z=0, where 
z is the proper length along the geodesies orthogonal to the wall. In order to obtain a thick 
gravitating wall solution one should solve the coupled system of the Einstein and scalar field 
equations which are as follows 
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(5) 



where is the Ricci tensor. For a static field X(z), the Einstein equations (@J) leads to the 
constraint B(z,t) = A(z) exp(fct). 

BCG investigated the spacetime of a thick gravitating planar domain wall for a A</> 4 potential 
|17j . In this context, the dimensionless parameter e in front of the equation @ can be taken to 
characterize the coupling of gravity to the scalar field: 
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(6) 



Supposing that gravity is weakly coupled to the scalar field, A(z) and X(z) may be expanded 
in the powers of e: 



A t (z) = A (z) + eA 1 (z) + O(e 2 ), 
X(z) = X (z) + eX 1 (z) + O(e 2 ), 



(7) 
(8) 



where we have introduced the index i to indicate the quantities inside the wall. 
In the limit e — > 0, one should obtain the same results as that of the non-gravitating planar wall, 
i.e. Ai(z) = 1 and X(z) = tanh(^). Using these expansions, BCG solved the coupled Einstein 
and scalar field equations to first order in e and obtained the following results: 
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(9) 

(10) 
(11) 



Eqs. ()9I11|) represent the perturbative solution to the spacetime of the thick wall described by 
Eq. In the presence of gravity, there is an event horizon at the proper distance z = \ from 
the core of the wall where Ai(z) goes to zero and the coordinate system used in © breaks down 
as discussed by BCG. 
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3 The thick wall formalism 



In this section we first make a short review of the thick wall formalism developed by KM in 
|19| . Then we apply it to the thick planar domain wall described by the metric jHJ). In KM 
formalism a thick wall is modelled with two boundaries Si and £2 dividing a spacetime M 
into three regions. Two regions M+ and M_ on either side of the wall and region Mo within 
the wall itself. Treating the two surface boundaries £1 and £2 separating the manifold .Mo 
from two distinct manifolds A4 + and .M_, respectively, as nonsingular timelike hypersurfaces, 
we do expect the intrinsic metric huv and extrinsic curvature tensor of (j=l,2) to be 
continuous across the corresponding hypersurfaces. These requirements named the Darmois 
conditions are formulated as 

[habkj=0 J = 1,2, (12) 
[Kabhj =0 i = 1,2, (13) 

where the square bracket denotes the jump of any quantity that is discontinuous across £j, Latin 
indices range over the intrinsic coordinates of Ej , and Greek indices over the coordinates of the 
4-manifolds. To apply the Darmois conditions on two surface boundaries of a given thick wall 
one needs to know the metrics in three distinct spacetimes M+, A4- and .Mo being jointed at 
While the metrics in M. + and Ai- are usually given in advance, knowing the metric whin 
the wall spacetime Aio requires a nontrivial work. 

Let us now impose these junction conditions for the self gravitating thick planar domain wall 
described in the previous section. Recalling e<l, corresponding to weak self gravity, we first 
introduce a dimensionless parameter 1 > A > 1 to assume the scalar field takes its vacuum 
values on the wall boundaries Si and £2 being located at the proper distances z = ± Aw/2 far 
from the wall's core surface while the coordinate system in (J3J) is still valid. We may then think 
of Aw as the proper thickness of the planar domain wall. Now, in the coordinate frame of the 
metric (j3J), in which the wall is stationary, the nonvanishing components of the intrinsic metric 
h^y and extrinsic curvature of £.,• take the following simple forms: 

h^u = g^u, V,v±z, (14) 

K iw = -\dtiu,z- (15) 

In order to find the spacetime metric on both sides of £j, we first note that within the vacuum 
region M+ (-M_) in which <f> = rj(—rj), the spacetime metric can be easily determined by solving 
the Einstein equations @ yielding 

A (z) = -k \z\ + C Q , (16) 

where the index o denotes the quantities in the vacuum region, and k Q and C Q are the relevant 
integration constants. Using Eqs. (|14j) and (|15|) we write down the junction conditions (|12|) and 
(fTS]) as 

k = k a , (17) 
M z )\z=Aw/2 = A (z)\ z=Aw/2 , (18) 
dMz) _ dA Q (z) 

dz \z=Aw/2 ~ g z \z=Aw/2- liyj 
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We now use the solutions ©, lfTU|l and (|TB|) due to BCG for the wall metric in the region 
Mo, where the scalar field varies according to (jllj) . and for the metric in the vacuum regions, 
respectively. Then the junction conditions (|18|) and (|19|) lead to the following constraints on the 
vacuum metric constants C Q and k Q 

Co =l + k ^- e -{2 ln(cosh + l tanh 2 ^ , (20) 

k Q = — ( tanh —-- tanh 3 — \ . (21) 
w \ 2 3 2 J v ' 

Note that within the context of BCG work it is supposed that the boundaries of the wall, where 
the scalar field takes its vacuum values, are at infinity. In contrast, we have modelled the thick 
planar wall in such a way that the wall boundaries Sj are situated at the finite proper distances 
±Aro/2 from the core of the wall. Hence, we choose A to be sufficiently large in order to 
simulate BCG solution within our thick wall model. 



4 From the thick to thin domain walls 

We now turn our attention to the thin wall limit of our thick wall model. First let us define the 
process of passing from a thick gravitating domain wall to a thin one by letting e and w go to 
zero while keeping their ratio ^ fixed. This has the effect that the distance of the event horizon 
to the domain wall remains finite. The continuity of the extrinsic curvature tensor across 
the thick wall boundary, say Si, located at the proper distance z = Au>/2 as a consequence of 
the Darmois junction conditions (|13|) yields: 

K° l/ \ z= Aw/2 = K^ u \ z=Aw / 2 . (22) 

Considering {xx) component of the condition (|22j) and using (|15|) we can evaluate the right hand 
side of Eq. (|22|) for the BCG wall metric solution ([§]). This yields 

K° x \ z=Aw /2 = ^ (tanh ^ - ^ tanh 3 ~ | ^21n(cosh —) + i tanh 2 exp(2k a t), 

(23) 

where we have also used the condition HI 7(1 . Imposing the above thin wall limit prescription, 
Eq. Q23JI reduces to 

K x \ z=0 = - (tanh ^ - 1 tanh 3 -) exp(2/c t), (24) 
w V 2 3 I J 

To identify the right hand side of the equation ((24(1 we recall the definition of the surface energy 
density a of an infinitely thin wall. Within our thick wall model it takes the form 

/•Aw/2 

a = lim / pdz, (25) 

(w-^O.e^O) J -Aw/2 

where p = p(z) is the energy density of the scalar field determined by the BCG scalar field 
solution I|11JI . Using Eq. (|25|) we obtained the following expression for a: 

1 e / , A 1 , o AN 
CT = ^( tanh 2-3 tanh Tj' (26) 



6 



where we used the definition of e given by ©• Comparing the results (|24|) and Q26|) one imme- 
diately obtains 

K° x \ z=0 = 2irGaexp(2k t). (27) 

Now, this is just the (xx) component of the Israel junction condition K° u = 2ixGah ilv for a 
planer thin domain wall placed at z = 0. The Israel thin wall approximation treats the wall as 
a singular hypersurface with the surface energy a separating the two plane symmetric vacuum 
spacetimes A4 + and A4- from each other. We see, therefore, within our thick wall formulation, 
in the thin wall limit the Darmois junction conditions for the extrinsic curvature tensor at the 
wall boundaries generate the well-known Israel jump condition. In the process of passing from 
a thick planar domain wall to the thin one all the information about the internal structure of 
the wall is squeezed in the parameter a characterizing the wall surface energy density according 
to ©. 



5 Evolution of the thick planar domain wall 

In the coordinate system (J3J) used to describe the thick wall, the wall is locally static and its core 
is situated at the fixed position z = 0. To study the evolution of the thick wall we will look for 
a reference frame in which the wall is moving. In order to do this, let us introduce a coordinate 
system (v, r, X, Y) in which the plane symmetric metric © with the solutions © and for 
A and k, respectively is transformed to the form 

ds 2 = f{z, t)dv 2 + 2H(z, t)dvdr - r 2 {dX 2 + dY 2 ), (28) 

by applying the following appropriate coordinate transformations: 

X = kx, 
Y = ky, 

r _ M*) c kt 
k ' 

v = v(z,t), (29) 

where v(z,t) is an unknown function to be determined. Identification of the metrics (J5J and 
(fgSJ) leads to 

x / dv \ 2 , Bv dr , 9 , N 

><*•*>(«) + 2H ^aiai = A ^ 

n f dv\ 2 „. .dvdr 

. dv dv . ( dv dr dv dr \ , . 



Solving the above equations we find 



^• t » = ■ (31) 
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and 



at h U2 ( * )2 _ ( & )a 



ft; A / Agg - g2 gj 



a * # U 2 (&) 2 -(&) 2 



(32) 



where the sign parameters s% and S2 can be taken independently to be ±1. Focusing our attention 
on the z > side of the wall, we infer that the continuity of the extrinsic curvature on the wall 
boundary placed at z = Aw/2 , however, requires si = — 1 and S2 = — 1. 
In addition, for the coordinate transformations to be integrable we require: 



d v d 2 v 



dtdz dzdt 

Imposing the condition (|33j) on the solutions (|32|) we arrive at the following: 

1 dH(z,t) A(z) dH(z,t) 



(33) 



H(z,t) dt H(z,t) dz 



9(z), (34) 



M=- A 'u- k+ W%- (35) 

We have now all the prerequisites to find the unknown functions defining the transformation 
and to understand the dynamics of the thick wall. Assuming a factorizable solution in the form 
H(z,t) = h z (z)ht(t), then Eq. (|34jl leads to the two following equations: 

1 ^ = C, (36) 



h t (t) dt 
A(z) dh z {z 



h z (z) dz 



+ 9(z) = C, (37) 



where C is an arbitrary constant. Integrating Eqs. (|36|) and (|37|) over t and z, respectively, we 
find 

H(z, t) = Dh z (z) exp(Ct), (38) 

with 

2 e (C+k)*-ftanh 2 £ 

hz(z) = 27, (39) 

(2-tanhf)(l+tanhf)2(coshf)- 

where D is an integration constant. The transformation (|29|) is now fully determined, except 
for the explicit functional dependence of v which will be considered later. The wall in the new 
coordinates is defined by R(t) = r(t(r), z = const) with r being its proper time as measured by 
observers at rest on the z = const surface. It is then easily seen that dynamics of the thick wall 
in the new coordinates is given by 

A(z) kT 

R(r) = -p e ^T . (40) 

K z=const 
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This dynamical equation takes a simple form for the core of the wall defined by R(tq) = 
r(t(ro), z = 0). In this case we have 

R(ro) = \e kT °. (41) 

Comparing Eq. (|41|) with the thin wall dynamics of [S], it is obvious that the core of our thick 
wall evolves as if it were a thin wall with the effective surface tension a = -^q- in addition, 
taking into account Eqs. H17|) . H21|) . and (|26[) it is seen that in the thin wall limit Eq. (|41j) is 
reduced to the Ipser-Sikivie solution for the evolution of a thin planar wall |Hj. 
We now intend to express the equation of motion (|41|) as a function of the new coordinate v. 
Note that from the 3-metric h a b induced on each z = const surface of the wall we have: 

f(z, r)i 2 (z, t) + 2R(t)H(z, t)v(z, t) = 1, (42) 

z=const 

where dot denotes the derivative with respect to r. Solving this for v gives 

(43) 



RH + JR 2 H 2 + / 



/ 



z=const 



where we have made the sign choice to make sure of having a finite solution on each z = const 
surface specifically on the wall boundaries where the continuity of the extrinsic curvature requires 
/ given by the expression lj31|) to be zero. To proceed further, we focus our attention on the 
core of the thick wall and first write down explicitly the solution (|l3*)) there 



e -(C+k)r 

D 



(44) 



where we have used Eqs. (|3*Tj) . (J3HJ), (|39j) . and (|4*T|) . Now, requiring that tq and v point to the 
same direction, i.e. v > 0, leads to the constraint D > 0. Integrating this equation over To we 
get the following results 

<r°) = O^{l-e- {C+k)T0 )> C>-k, (45) 

^o) = n ( e ' c+fc|ro " A . c < ( 46 ) 



D\C + k\ 

where we have put v = at tq = 0. Inserting the solutions ()45|) and (|46|) into Eq. ()41|) yields 
the functional dependence of the equation of motion of the wall's core on the time coordinate v 
as follows: 

R(v) =-(l-D(C + k)v)v& , O-k, (47) 

Ki 

R(v) = - (1 + D\C + k\v)W+® , C<-k. (48) 
k 

It follows that the core of the wall expands when looked upon from the bulk frame v regardless 
of whether C > k or C < k. 
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To learn more about the evolution of the wall, we write down the derivatives of Eqs. (|47|) and 
(148(1 with respect to the coordinate time v 



2k + C 

^ = D f i n,n K ) ™ > C > ~ k > (49) 

— = D (1 + D\C + k\v)\^\~\ C<-k. (50) 
ai> 

Consequently, from Eqs. (|49|) and (|50|) we see that the expansion of the wall's core slows down 
provided C < —2k, otherwise it speeds up. 

From the coordinate transformation we have chosen to investigate the dynamics of the wall, it 
is clear that the wall thickness is no longer a constant . To see the time dependence of the 
thickness we first note that the spacetime parts (z > 0) and (z < 0) are transformed separately 
to spheres. Therefore, it is more appropriate to look at the half-thickness (w*) defined by 

w* =R(v)\ z=0 -R(v)\ z= A m . (51) 

Following the above procedure leading to the equations of motion for the core of the wall one 
can find the corresponding equation for the boundary of the wall located at z = We get 

. . . At 2D(C + k)h z (z) , A , 
R{v)\ z= ^ = -[\ 1 A >z{) (v-vo)) L =¥ , (52) 

where we have considered the case C > —k, vo denotes the initial value of v when the proper 
time on the boundary of the wall is taken to be zero. Plugging Eqs. (|47|) and (|52|) into the 
definition (|51|) we end up with 

w* = - (1 - D(C + k)v)v& 



k 



-k 



At 2D{C + k)h z (z) ( A^, 

T 1 - A (v-vq)) | ah, (53) 



k V A V ' z= — 

where A{z) and h z (z) are given by the solutions © and (|39|) . respectively. This gives the time 

dr 
(ir. 



dependence of the half-width of the wall in terms of the new coordinate v. Noting that §f- < 0, 



we can see that the half-width w* is an increasing function of v. 



6 Conclusion 

We have studied the thick planar domain wall as a spacetime having two boundaries at the 
same proper distance from the wall's core, as first formulated in |19|. It has then been shown 
that in the thin wall limit the Darmois junction conditions generate the familiar Israel's jump 
condition at the boundary of the corresponding thin wall with the embedding spacetimes. It is 
realized that in the limiting process to the thick planar domain wall all the information about 
the internal structure of the wall is squeezed in the surface energy density a introduced in (|2(i|) . 
We have also given explicitly a coordinate transformation from the plane symmetric metric 
describing a thick planar domain wall at rest to a reference frame in which the thick wall evolves. 
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Starting from an ansatz for the transformed metric of the form (|28|) this aim was achieved by 
analytically solving the transformation equations. It is then seen that the core of the thick wall 
is moving with respect to the transformed reference frame as if it were a thin wall. In contrast, 
according to Eq. (|4U|) the other layers of the wall evolve differently. We have also studied the 
thin wall limit of the moving thick wall and showed that in this limit the thick wall's solution 
becomes the well-known thin wall solution of [SJ. In this dynamical picture, the wall thickness 
shows a time variation, which has been given explicitly. 
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